Evolution of the calcium feedback steps of vertebrate phototransduction

We examined the genes encoding the proteins that mediate the Ca-feedback regulatory system in vertebrate rod and cone phototransduction. These proteins comprise four families: recoverin/visinin, the guanylyl cyclase activating proteins (GCAPs), the guanylyl cyclases (GCs) and the sodium/calcium-potassium exchangers (NCKXs). We identified a paralogon containing at least 36 phototransduction genes from at least fourteen families, including all four of the families involved in the Ca-feedback loop (recoverin/visinin, GCAPs, GCs and NCKXs). By combining analyses of gene synteny with analyses of the molecular phylogeny for each of these four families of genes for Ca-feedback regulation, we have established the likely pattern of gene duplications and losses underlying the expansion of isoforms, both before and during the two rounds of whole-genome duplication (2R WGD) that occurred in early vertebrate evolution. Furthermore, by combining our results with earlier evidence on the timing of duplication of the visual G-protein receptor kinase genes, we propose that specialization of proto-vertebrate photoreceptor cells for operation at high and low light intensities preceded the emergence of rhodopsin, which occurred during 2R WGD.This work was supported by the Australian Research Council (grant nos. CE0561903 and DP110103294)

Quantitative modeling of the molecular steps underlying shut-off of rhodopsin activity in rod phototransduction

PURPOSE: To examine the predictions of alternative models for the stochastic shut-off of activated rhodopsin (R*) and their implications for the interpretation of experimentally recorded single-photon responses (SPRs) in mammalian rods. THEORY: We analyze the transitions that an activated R* molecule undergoes as a result of successive phosphorylation steps and arrestin binding. We consider certain simplifying cases for the relative magnitudes of the reaction rate constants and derive the probability distributions for the time to arrestin binding. In addition to the conventional model in which R* catalytic activity declines in a graded manner with successive phosphorylations, we analyze two cases in which the activity is assumed to occur not via multiple small steps upon each phosphorylation but via a single large step. We refer to these latter two cases as the binary R* shut-off and three-state R* shut-off models. METHODS: We simulate R*’s stochastic reactions numerically for the three models. In the simplifying cases for the ratio of rate constants in the binary and three-state models, we show that the probability distribution of the time to arrestin binding is accurately predicted. To simulate SPRs, we then integrate the differential equations for the downstream reactions using a standard model of the rod outer segment that includes longitudinal diffusion of cGMP and Ca²⁺. RESULTS: Our simulations of SPRs in the conventional model of graded shut-off of R* conform closely to the simulations in a recent study. However, the gain factor required to account for the observed mean SPR amplitude is higher than can be accounted for from biochemical experiments. In addition, a substantial minority of the simulated SPRs exhibit features that have not been reported in published experiments. Our simulations of SPRs using the model of binary R* shut-off appear to conform closely to experimental results for wild type (WT) mouse rods, and the required gain factor conforms to biochemical expectations. However, for the arrestin knockout (Arr−/−) phenotype, the predictions deviated from experimental findings and led us to invoke a low-activity state that R* enters before arrestin binding. Our simulations of this three-state R* shut-off model are very similar to those of the binary model in the WT case but are preferred because they appear to accurately predict the mean SPRs for four mutant phenotypes, Arr+/−, Arr−/−, GRK1+/−, and GRK1−/−, in addition to the WT phenotype. When we additionally treated the formation and shut-off of activated phosphodiesterase (E*) as stochastic, the simulated SPRs appeared even more similar to real SPRs, and there was very little change in the ensemble mean and standard deviation or in the amplitude distribution. CONCLUSIONS: We conclude that the conventional model of graded reduction in R* activity through successive phosphorylation steps appears to be inconsistent with experimental results. Instead, we find that two variants of a model in which R* activity initially remains high and then declines abruptly after several phosphorylation steps appears capable of providing a better description of experimentally measured SPRs.This work was supported by award number R01EY023603 from the US National Eye Institute

Analysis of Paralogons, Origin of the Vertebrate Karyotype, and Ancient Chromosomes Retained in Extant Species

A manually curated set of ohnolog families has been assembled, for seven species of bony vertebrates, that includes 255 four-member families and 631 three-member families, encompassing over 2,900 ohnologs. Across species, the patterns of chromosomes upon which the ohnologs reside fall into 17 distinct categories. These 17 paralogons reflect the 17 ancestral chromosomes that existed in our chordate ancestor immediately prior to the two rounds of whole-genome duplication (2R-WGD) that occurred around 600 Ma. Within each paralogon, it has now been possible to assign those pairs of ohnologs that diverged from each other at the first round of duplication, through analysis of the molecular phylogeny of four-member families. Comparison with another recent analysis has identified four apparently incorrect assignments of pairings following 2R, along with several omissions, in that study. By comparison of the patterns between paralogons, it has also been possible to identify nine chromosomal fusions that occurred between 1R and 2R, and three chromosomal fusions that occurred after 2R, that generated an ancestral bony-vertebrate karyotype comprising 47 chromosomes. At least 27 of those ancestral bony-vertebrate chromosomes can, in some extant species, be shown not to have undergone any fusion or fission events. Such chromosomes are here termed "archeochromosomes," and have each survived essentially unchanged in their content of genes for some 400 Myr. Their utility lies in their potential for tracking the various fusion and fission events that have occurred in different lineages throughout the expansion of bony vertebrates

Long-term storage and age‐biased export of fluvial organic carbon: field evidence from West Iceland

Terrestrial organic carbon (OC) plays an important role in the carbon cycle, but questions remain regarding the controls and timescale(s) over which atmospheric CO₂ remains sequestered as particulate OC (POC). Motivated by observations that terrestrial POC is physically stored within soils and other shallow sedimentary deposits, we examined the role that sediment storage plays in the terrestrial OC cycle. Specifically, we tested the hypothesis that sediment storage impacts the age of terrestrial POC. We focused on the Efri Haukadalsá River catchment in Iceland as it lacks ancient sedimentary bedrock that would otherwise bias radiocarbon‐based determinations of POC storage duration by supplying pre‐aged “petrogenic” POC. Our radiocarbon measurements of riverine suspended sediments and deposits implicated millennial‐scale storage times. Comparison between the sample types (suspended and deposits) suggested an age offset between transported (suspended sediments) and stored (deposits) POC at the time of sampling, which is predicted by theory for the sediment age distribution in floodplains. We also observed that POC in suspended sediments is younger than the predicted mean storage duration generated from independent geomorphological data, which suggested an additional role for OC cycling. Consistent with this, we observed interparticle heterogeneity in the composition of POC by imaging our samples at the microscale using X‐ray absorption spectroscopy. Specifically, we found that particles within individual samples differed in their sulfur oxidation state, which is indicative of multiple origins and/or diagenetic histories. Altogether, our results support recent coupled sediment storage and OC cycling models and indicate that the physical drivers of sediment storage are important factors controlling the cadence of carbon cycling

Human scotopic dark adaptation: Comparison of recoveries of psychophysical threshold and ERG b-wave sensitivity

We have compared the time course of dark adaptation of the human scotopic visual system, measured psychophysically and from the b-wave of the electroretinogram (ERG), for bleaches ranging from a few percent to near total. We also measured light adaptation, in order to apply a "Crawford transformation" to convert the raw measurements of dark adaptation into equivalent background intensities. For both the "psychophysical threshold equivalent" intensity and the "ERG b-wave sensitivity equivalent" intensity, the equivalent background declined over much of its range with an "S2" component, though with somewhat different slopes of -0.36 (psychophysical) and -0.22 (ERG) log10 unit min-1, respectively. In addition, the magnitude of the equivalent background was approximately 1 log10 unit lower in the psychophysical experiments than in the ERG experiments. Despite these differences, the two approaches extract a common time course for the decline in level of free opsin following moderately large bleaches. We conclude that the recovery of psychophysical scotopic visual threshold over the S2 region reflects events that are present by the stage of the first synapse of rod vision, stemming ultimately from the presence of unregenerated opsin in the rod outer segments

Evolution of the calcium feedback steps of vertebrate phototransduction

We examined the genes encoding the proteins that mediate the Ca-feedback regulatory system in vertebrate rod and cone phototransduction. These proteins comprise four families: recoverin/visinin, the guanylyl cyclase activating proteins (GCAPs), the guanylyl cyclases (GCs) and the sodium/calcium-potassium exchangers (NCKXs). We identified a paralogon containing at least 36 phototransduction genes from at least fourteen families, including all four of the families involved in the Ca-feedback loop (recoverin/visinin, GCAPs, GCs and NCKXs). By combining analyses of gene synteny with analyses of the molecular phylogeny for each of these four families of genes for Ca-feedback regulation, we have established the likely pattern of gene duplications and losses underlying the expansion of isoforms, both before and during the two rounds of whole-genome duplication (2R WGD) that occurred in early vertebrate evolution. Furthermore, by combining our results with earlier evidence on the timing of duplication of the visual G-protein receptor kinase genes, we propose that specialization of proto-vertebrate photoreceptor cells for operation at high and low light intensities preceded the emergence of rhodopsin, which occurred during 2R WGD

Evolution of the vertebrate phototransduction cascade activation steps

We examine the molecular phylogeny of the proteins underlying the activation steps of vertebrate phototransduction, for both agnathan and jawed vertebrate taxa. We expand the number of taxa analysed and we update the alignment and tree building methodology from a previous analysis. For each of the four primary components (the G-protein transducin alpha subunit, GαT, the cyclic GMP phosphodiesterase, PDE6, and the alpha and beta subunits of the cGMP-gated ion channel, CNGC), the phylogenies appear consistent with expansion from an ancestral proto-vertebrate cascade during two rounds of whole-genome duplication followed by divergence of the agnathan and jawed vertebrate lineages. In each case, we consider possible scenarios for the underlying gene duplications and losses, and we apply relevant constraints to the tree construction. From tests of the topology of the resulting trees, we obtain a scenario for the expansion of each component during 2R that accurately fits the observations. Similar analysis of the visual opsins indicates that the only expansion to have occurred during 2R was the formation of Rh1 and Rh2. Finally, we propose a hypothetical scenario for the conversion of an ancestral chordate cascade into the proto-vertebrate phototransduction cascade, prior to whole-genome duplication. Together, our models provide a plausible account for the origin and expansion of the vertebrate phototransduction cascadeThis work was supported by the Australian Research Council (Grants CE0561903 and DP110103294). Sequences have been deposited in GenBank with accession numbers KY820586 – KY820625

A quantitative account of mammalian rod phototransduction with PDE6 dimeric activation: Responses to bright flashes

We develop an improved quantitative model of mammalian rod phototransduction, and we apply it to the prediction of responses to bright flashes of light. We take account of the recently characterized dimeric nature of PDE6 activation, where the configuration of primary importance has two transducin molecules bound. We simulate the stochastic nature of the activation and shut-off reactions to generate the predicted kinetics of the active molecular species on the disc membrane surfaces, and then we integrate the differential equations for the downstream cytoplasmic reactions to obtain the predicted electrical responses. The simulated responses recover the qualitative form of bright-flash response families recorded from mammalian rod photoreceptors. Furthermore, they provide an accurate description of the relationship between the time spent in saturation and flash intensity, predicting the transition between first and second ‘dominant time constants’ to occur at an intensity around 5000 isomerizations per flash, when the rate of transducin activation is taken to be 1250 transducins s−1 per activated rhodopsin. This rate is consistent with estimates from light-scattering experiments, but is around fourfold higher than has typically been assumed in other studies. We conclude that our model and parameters provide a compelling description of rod photoreceptor bright-flash responses.This work was supported by award no. R01EY023603 from the US National Eye Institute

Evolution of vertebrate retinal photoreception

Recent findings shed light on the steps underlying the evolution of vertebrate photoreceptors and retina. Vertebrate ciliary photoreceptors are not as wholly distinct from invertebrate rhabdomeric photoreceptors as is sometimes thought. Recent information on the phylogenies of ciliary and rhabdomeric opsins has helped in constructing the likely routes followed during evolution. Clues to the factors that led the early vertebrate retina to become invaginated can be obtained by combining recent knowledge about the origin of the pathway for dark re-isomerization of retinoids with knowledge of the inability of ciliary opsins to undergo photoreversal, along with consideration of the constraints imposed under the very low light levels in the deep ocean. Investigation of the origin of cell classes in the vertebrate retina provides support for the notion that cones, rods and bipolar cells all originated from a primordial ciliary photoreceptor, whereas ganglion cells, amacrine cells and horizontal cells all originated from rhabdomeric photoreceptors. Knowledge of the molecular differences between cones and rods, together with knowledge of the scotopic signalling pathway, provides an understanding of the evolution of rods and of the rods' retinal circuitry. Accordingly, it has been possible to propose a plausible scenario for the sequence of evolutionary steps that led to the emergence of vertebrate photoreceptors and retina

Implications of dimeric activation of PDE6 for rod phototransduction

We examine the implications of a recent report providing evidence that two transducins must bind to the rod phosphodiesterase to elicit significant hydrolytic activity. To predict the rod photoreceptor's electrical response, we use numerical simulation of the two-dimensional diffusional contact of interacting molecules at the surface of the disc membrane, and then we use the simulated PDE activity as the driving function for the downstream reaction cascade. The results account for a number of aspects of rod phototransduction that have previously been puzzling. For example, they explain the existence of a greater initial delay in rods than in cones. Furthermore, our analysis suggests that the ‘continuous’ noise recorded in rods in darkness is likely to arise from spontaneous activation of individual molecules of PDE at a rate of a few tens per second per rod, probably as a consequence of spontaneous activation of transducins at a rate of thousands per second per rod. Hence, the dimeric activation of PDE in rods provides immunity against spontaneous transducin activation, thereby reducing the continuous noise. Our analysis also provides a coherent quantitative explanation of the amplification underlying the single photon response. Overall, numerical analysis of the dimeric activation of PDE places rod phototransduction in a new light.This work was supported by award no. R01EY023603 fromthe US National Eye Institut